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1. Introduction 

Quantum field theories have come a long way since its inception and today, all known 
interactions have been classed within its model description for elementary particles in 
nature Within this framework, studies on Feynman loop integrals became even more CIT: giover-qcd 
compelling, challenging us with ever increasing mathematical complexities associated with 
the perturbative approach. Results for calculations of diagrams with massive internal 
particles were presented by the authors j|] and others J|]; non-planar double-box scalar cm equiv 
and tensorial integrals were studied by Tausk [J|] and Smirnov using Mellin-Barnes CIT: tramontano 
technique — recently the triple box as well — while Gehrmann and Remiddi || using ^ 
the powerful differential equation method calculated several 2-loop integrals. Bern and q T 
collaborators studied dimensionally regularized one-loop pentagon integrals Q as well as ciT: bem-penta 
Binoth et al [^] (numerically) and the authors (analytically and with arbitrary exponents CIT: binoth 
of propagators) tackle one-loop scalar hexagon integrals || . cm n-point 

These three methods are very powerful and interesting: integration-by-parts relates a 
complicated integral to simpler ones with some exponents of propagators raised to powers 
greater than one. Differential equation method also relates a complicated graph to simpler 
ones - that means lesser number of loops or simpler graphs with the same number of loops -, 
with the exponents of propagators also changing when one uses this method. Mellin-Barnes 
approach is based on an integral representation in which the resulting integrals can be 
carried out summing up the residues (either in the right or left complex plane) . The latter 
one therefore yields the two domains of validity available for the hyper geometric functions 
connected by analytic continuation. These in turn define the two distinct kinematical 
regions of interest [10]. cm boos 



An altogether different and elegant approach has been suggested back in the middle 



of the 1980's by Halliday and Ricotta [11], coined NDIM, where the key point is the cm haiiiday 
introduction of a negative dimensional space to work out the integral. Their seminal ideia 
was reframed within the context of solving systems of linear algebraic equations and from 
their original work to our present understanding and experience in it we deem NDIM to be 
a more compleat tool to handle complex Feynman integrals, covariant and non-covariant 
alike, and far simpler in its essence to implement, needing only to deal with systems of linear 
algebraic equations of first degree. All kinematical regions of interest come simultaneously 
defined and the bonus by-product is that it defines even as yet unknown relationships of 
analytic continuation among hyper geometric functions of several variables. 

Our aim in this work is to evaluate some loop integrals which were not considered in the 
literature up to now using the NDIM technique. They are double-box scalar integrals with 
six propagators: in the simplest case, i.e., massless internal particles and on-shell external 
legs the result j|] is well-known for a special case of exponents of propagators (all equal cm tausk 
to minus one); we fill the gap presenting the result for arbitrary exponents of propagators 
and extend our knowledge studying the same graph where four internal particles have 
arbitrary mass \i and also studying the case where one of the external legs is off-shell. 
These diagram computations become important as progress in perturbative calculations 
for fundamental interactions between particles are checked against the background of our 
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Figure 1: Scalar non-planar double-box with four massive propagators which are represented by 
thick lines. The labels in the internal lines represent the exponents of propagators. All external 
momenta are considered to be incoming, so pi + P2 + P3 + Pa = 0. 



present experimental data increases their precision measurements. 

The outline for our paper is as follows: in section 2 we study covariant four-point 
integrals with six propagators - non-planar double-box with four massive propagators - the 
exact result of which is written in terms of hyper geometric function of three variables. We 
also perform numerical calculations in order to expand the result in powers of e. Section 3 
is concerned with massless non-planar double-box with off-shell external legs and in section 
4 we give concluding remarks for our present work. 



2. Massive non-planar double-box 

To begin with consider the integral for the massive non-planar double-box with six propa- 
gators, namely, 



'-M 



d D q d D r ( (? 2 V) i [( g -p 3 )^] J [(9+r)'V] fe [(9+r+P2)'V] < (^) m (r-p4) 2n • (2.1) LAB: i 



\2 ,,2-ifc 



2 ,.2-iJ/ „2\mt 



This is the scalar integral which arises in the computation of the diagram depicted in 
figure 1, where the external legs represent on-shell massless particles, i.e., p\ = p\ = P 2 = 
pj = 0. 

The generating functional for massless on-shell non-planar double-box ll^] is given by, CIT: pentabox 



Go = J d D qd D r exp {-aq 2 - f3(q - p 3 ) 2 --f(q+r) 2 -8(q + r+p 2 ) 



2 i 2 

<pr 



u(r -Pi) 2 } 

TT 2 \ D/2 

- exp 



(P-ftus + adujt + (30(f>u) 



(2.2) 

(2.3) LAB: geradora-on 
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where we have defined A' = a + (3 + 7 + 9, A = 07 + aO + (3j + (39 + \'{<p + to) and s, t, u 
are the usual Mandelstam's variables, given by 

s = 2p 1 -p 2 , t = 2p 1 -p 3 , u = 2p 1 -p 4 . (2.4) 

Since pf = 0, (i = 1, 2, 3, 4), observe that s + t + u = follows from the above equation. 
Then, considering some of the internal particles to be massive, that is, diagram of 
figure 1, one has the generating functional, 

G M = exp(\'n 2 )G , (2.5) 

with the massive sector factorized, and the following system of algebraic equations, 

X 2 + Y1234 + Wi = i 

X13 + Y 5678 + W 2 = j 

Xx + Y159 + W 3 + Z 1 = k 

X 2 3+Y 2 e + W 4 + Z 23 = I , (2.6) LAB:sistema 

X3 + Y 379 + Z 2 = m 
X12 + I48 + Z13 = n 

sx + sy + sz = -D/2 

where Wj are the indices labelling the pure massive sector. We use the shorthand notation, 

X a bc = X a + Xfr + X c , 

and so on. In the last equation EX = X i23 , YY = ^1234567895 = ^123- Note that 
the total of sum indices are 19 with 7 constraint equations, so that the result will be 
a series of 12 indices. This 12-fold sum may be constructed in various ways, in fact, 
C 7 12 ways. A huge number of ways. The majority of them - 30,972 - yield vanishing 
determinant for the algebraic system, so that the solution is trivial in all of these cases. 
The remaining 29,416 yield non-vanishing determinant and non-trivial solutions expressed 
as hyper geometric series representations. The several variables that identify these series 
belong to a subset of the following set 

tusstut u s Ap? 4fi 2 4fi 2 "1 
' ? I' «' ? «' ? V' 4^2' 4^2' ~T' V T" J ' { ' 

The simplest of the hypergeometric series representations for 2m is given by a triple 
series, 

T = D ( 2 Y V V X^Y^Z^(-a\X l2 3){-n\X l2 )(-m\X3){-l\X 2 3) 
M x u ^tx 3 =o X 1 !X 2 !X 3 !( J D/2|X 123 )(-i-j|X 123 ) 

„ (-ilXzX-jlXuX-klXM-k-l-m-n-D/qXm) 

{Z.o ) LAB: tnpla 



(1/2 - a/2 - m/2 - n/2|X 123 )(-a/2 - m/2 - n/2|X 123 ) ' 

where a = i+ j + k + l + m + n + D, is the sum of exponents and dimension, (a\b) is the 
Pochhammer symbol, 

/ i7 \ / \ T(a + 6) 

(a\b) = (a) b = -^- r , 
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and 



(— a — m — n\m + n){D /2\m + n){— i — j\ — m — n — D/2)(—k — l\ — m — n — D/2), (2.9) 



where the subset of three variables are 



X 



s 



Y 



t 



z 



u 



Observe that the final result, eq. (2.8), has only a three-fold series. However, the 
expression provided by the solution of system ( |2.6D was a 12-fold series. It is very easy 
to understand why this is so. Among the defining variables for the hypergeometric series 
representations here there are 9 whose variables are just unity and are summable series. 



In other words, we were able to sum up nine of them using Gauss' summation formula[13]. CIT: luke 
Our strategy is therefore to choose as many series as possible in which the individual sums 
can be written as 2 iq (a, b; c|l), then we plug them in a computer program that do the 
job, i.e., sum up the series using Gauss' summation formula. 

Once we have evaluated this result, it is a straightforward exercise to write down several 
other solutions which are connected by symmetry in the diagram, namely, by exchanging 
the pairs (i «-> k, s «-► t), and (j «-► I, s «-► t) 

If we are interested in the primary integral where all the exponents of propagators are 
equal to minus one, the above result reduces to 



-M 



^V)^r({-i}) £ 

(1|X 1 )(1|X 13 )(1|X 2 ) 



x x iy x 2Z x 3(6 _ D|x 123 )(l|X 12 )(l|X 3 )(l|X 2 3) 



X 1 !X 2 !X 3 !( J D/2|X 12 3)(2|X 123 ) 



(9/2-D/2\X : 



123 J 



where 



r({-i}) 



T(D/2 - 2)r(D/2)r(6 - £>)T 2 (4 - D/2) 



(2.10) 



(2.11) 



LAB: 4mass 



T(8 - D) 

The original system of linear equations defines a 7 x 19 rectangular matrix. From it we 
can draw 50, 388 square submatrices of dimension 7 x 7, of which 30, 972 yield vanishing 
determinant, as already said before. Of course, among the 29,416 solvable solutions that 
remain, NDIM provides other kind of series, such as 5-fold and 7-fold series, i.e., hyper- 
geometric series representations with five and seven variables respectively (meaning seven 
and five summable series expressed as 2 -Fj.(a, 6; c\l) respectively). And all of them have 
symmetries among s, t and u, namely, 



(p2 *-* P4, j <-> I, i <-> k, s <-> u), 



(p 2 <-> ps, I n, k <-> m, t <-> s) 



(2.12) LAB: simetrias 



so, for each hypergeometric series representations provided by NDIM there are other two, 
also originated from the system of algebraic equations, which represent the same integral 
and can be transformed in the first using (|2T^ ). This is the case of (|2[ 
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Figure 2: Scalar massless non-planar double-box with six propagators. 



2.1 Numerical calculation 

Expansion in e for the integral fl2,l| ) can be obtained numerically to all orders since our re- 
sult ( |2.10|) is exact. Hypergeometric series converge very fast and eq. (|2.10|) can be truncated 



after just few terms. We consider two examples, namely, s 
and s = 2,t = —1, u = —1, [i = 4. 



-l,t 



-2,8 



-3,V 



25 



Table 



N 


a 


b 


c 


1 


-0.15877787089947 


0.1144504436314 


-0.4375363472289 


2 


-0.15905593258958 


0.1141346540498 


-0.4382967223214 


3 


-0.15904157598309 


0.1141585957166 


-0.4382477388046 


4 


-0.15904248795073 


0.1141566914709 


-0.4382515309963 


5 


-0.15904248795073 


0.1141568521924 


-0.4382512066856 


6 


-0.15904242717753 


0.1141568379462 


-0.4382512361932 


7 


-0.15904242673854 


0.1141568392587 


-0.4382512333925 


8 


-0.15904242677727 


0.1141568391341 


-0.4382512336666 


9 


-0.15904242677373 


0.1141568391462 


-0.4382512336391 


10 


-0.15904242677406 


0.1141568391450 


-0.4382512336419 


20 


-0.15904242677403 


0.1141568391451 


-0.4382512336417 



N=Number of terms of each series and coefficients of e-expansion for (2^1). Case-I: s = 
—l,t = —2,u = — 3,4/i 2 = 25. It is outside the physical region and therefore represents a 
numerical sample. The result is in the form - + b + ce, where a, b, c are given in the table. 
It is possible to obtain the e-expansion to all orders, since our result is exact. 
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Table 



N 


a 


b 


c 


1 


-0.16666202713362544092 


0.11111642984629631817 


-0.4557060127212616592 


2 


-0.16668058021298683828 


0.11109516323349484374 


-0.4557569968657018264 


3 


-0. 16668072454737026304 


0.11109492182787342273 


-0.4557574905733665544 


4 


-0. 16668072849754143706 


0.11109491356240756367 


-0.4557575070366868277 


5 


-0. 16668072856344380806 


0.11109491340169955725 


-0.4557575073610663256 


6 


-0. 16668072856499855487 


0.11109491339744859266 


-0.4557575073698739456 


7 


-0.16668072856503197638 


0.11109491339734860742 


-0.4557575073700873608 


8 


-0. 16668072856503277706 


0. 1 1 109491339734602934 


-0.4557575073700930307 


9 


-0.16668072856503279614 


0.11109491339734596400 


-0.4557575073700931786 


10 


-0.16668072856503279661 


0.11109491339734596228 


-0.4557575073700931826 


20 


-0. 16668072856503279662 


0.11109491339734596224 


-0.4557575073700931827 



N=Number of terms of each series and coefficients of e-expansion for (|2.1| ). Case-II: s = 
2,t = —l,u = — 1, jU = 4. Observe that hyper geometric series converges very fast: twenty 
terms for each of the three series provides 18 figures precision. The result is in the form 
- + b + ce, where a, b, c are given in the table. 



3. Massless Double-box 

Our method can also be used to study Feynman diagrams where external legs are off- 
shell[[L4|]. Consider for instance the massless non-planar double-box of figure 1. Let pf ^ cm probing 
0, (i = 1,2,3,4), that is, all external legs off-shell. The generating functional is more 
complicated, namely, 

Ga-off = G exp (-(39ujpl -a 2 p\- a 3 p| - a 4 p|), (3.1) 

where Go is the generating functional ( |2.3j ) for the on-shell massless diagram and 

a 2 = a^0 + a0(j> + (3^9 + ^9(j) + ^9uj (3.2) 
03 = a/?7 + af39 + a(3cp + otfiuj + j3^(j) 
a4 = a^/tj + a(j)u> + 139(f) + 7^0; + Qcjxo. 

In the previous case we used the result s + t + u = 0, however, for the case at hand 
Mandelstam variables are, 

s = {Pi + P2) 2 = (P3 + Pa) 2 (3.3) 

t = (p x +p 3 ) 2 = (P2+Pi) 2 
U = {p\+ Pif = {P2 + P-3,) 2 , 

and s + t + u = pi + p\ + p\ + p\ / 0. 

In the next subsections we will study some particular cases. See Table 
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Table: Number of systems, solutions and type of results 



Diagram 


4 Equal Masses 


Massless (on) 


Massless (I) 


Massless (II) 


System 


7 x 19 


7 x 15 


7 x 16 


7 x 20 


Total number 


50,388 


5,040 


11,440 


77,520 


Solutions 


29,416 


2,916 


4,632 


34,994 


Result 


Triple Series 


Double Series 


p( 3 ) 
A 


9-fold series 



3.1 One leg off-shell, p\ ^ 

There are two distinct cases to consider when one external leg is off-shell: p\ ^ and 
p| ^ 0, for j = 2, 3, 4, since the diagram is symmetric under the change 



P2 <-> P3, 



P2 <-> Pi, 



Observe that the vertex where p\ is attached is quartic and all the others are triple and 
can be interchanged leaving the diagram unchanged. 
Mandelstam variables must be rewritten as, 



s = p l + 2pi ■ p 2 
t = p\ + 2px ■ p 3 
u = p\ + 2pi • p 4 , 



(3.4) 



and s + t + u = p\ = M\. So the generating functional becomes slightly different, 



OFF 



2\ D/2 



cxp 



A 



(3.5) LAB: geradora-off 



Go exp 



A 



comparing (3.5) with fl2.3| ) we conclude that the original system gained only one variable, 
i.e., the former system (double-box with 4 legs on-shell) was 7x 15 and the present (double- 
box with 1 leg off-shell) is 7 x 16. The total number of solutions is now 11,440 being 6,£ 
trivial systems and we must deal with 4,632 possible ones. 

The simplest hypergeometric series representations for I\_ fj are triple series, 



l-o// 



P 1 Xl P 2 X2 P? [i (-a\X 123 )(-k\X 1 )(-i\X 2 )(-m\X 3 ) 
X!!X 2 !X 3 !(1 + I - a\Xi)(l + j - a\X 2 )(l + n - a\X 3 ) ' 



(3.6) LAB: tripla-loff-pl 



^1,2,3=0 

where Pi = s/Mf, P 2 = t/M^, P 3 = u/Mf, and following the usual approach for massless 
diagrams in the NDIM context we have summed up 6 series. The above triple series is a 
Lauricella's function of three variables [jisfl, namely, 



CIT: luke 



l-off 



n D hF 



(3) 



—a; —k, —i, —m 
1 + / — a, 1 + j — a, 1 + n — a 



Pl,P2,P 3 



(3.7) 
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where for convergence |Pi| + |P 2 | + I-P3I < 1- We also define, 

ft = (Ml) a ((T + D/2\ - 2a - D/2)(i + j + m + n + D\ -m-n-D/2) 

x(k + l + m + n + D\ - k - I - D/2)(i + j + k + I + D\ - i- j-D/2), 

x (-ik)HH(-™H, 

the symmetries of the diagram 

i <-> k, j <-> I, s <-> t; i <-> m, j «-> n, s «-> u; «-> m, Z <-> n, t <-> u, 

are expressed in our final result. 
3.1.1 Special Case 

In the special case where i = j = k = l = m = n= — lwe have, 



(3.8) 



(3.9) 



X 1 - o// = vr D / 1 ({-l})^ 



(3) 



6 -D; 1,1,1 
6-D,6-A6-D 



with 



/i({-l}) = (Pi 



2 , 6 _ D r 3 (L>- 5)r 3 (D/2 -2)r(6-D) 



(3.10) 



(3.11) 



r 3 (z?-4)r(3Z?/2-6) 

which have a double pole in the D = 4 — 2e limit. 

If one were interested in writing the above result in terms of more complicated func- 
tions, namely, logarithms, polylogarithms and S a> b integrals, then the following integral 
representation for is in order 



,(3) 



L 7, 7', 7" 



Pl,P2,P 3 

x(i- X2 y'-P'- 1 (i-x 3 ) 

where the integration is constrained to x\ + x 2 + X3 = 1 and 



! s /J ~V" 1 a^"" 1 ri-a;i^-' 9 - 1 



(1 - xiPi - x 2 P 2 - x 3 P 3 ) c 

7 "-/3"-l 



(3.12) 



r(/3)r(/?')r(/3")r( 7 - W{i - PW - p") 



in the present case one get, 



F 



(3) 



6 -D; 1,1,1 
6 - D,6 - D,6-D 



Pi,P2,P 3 



r( 7 )r(y)r(Y' 



— / dxidx 2 dx ? ,- r 

Ta Jo I 1 - a^i-Pi - x 2?2 - x 3 P 3 ) b u 



x(l-x 2 ) 4 -^(l-x 3 ) 



\4-D 



(3.13) 



then one take D = 4 — 2e and uses Taylor expansion. However, to carry out the integral 
of second derivatives of such integral representation can not be an easy task. For this 
reason we claim hypergeometric series representations are simpler than the ones in terms of 
polylogarithms: fast convergence, compact expressions and analytic continuation relations 
among (i.e. kinematical regions) them. 
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3.2 One leg off-shell, p| ^ 



Now we turn to the last case, the one where the external leg attached to a quartic vertex 
is off-shell. The generating functional is, 



Gi-off = G exp 



crvfl + a6(j) + 876 + 766 + j6lo 9 
A P > 



(3.14) LAB: geradora-p2 



we see immediately that there will be four extra sums. Simply compare (3.14) and (|3.5|). 
the former has four arguments (which properly expanded in Taylor series will produce the 
referred extra sums) more than the latter. 

The system of algebraic equations will be slightly different than fl2.6| ), i.e., for the 
present case one have, 



X2 + Y1234 + U12 

X13 + ^5678 + U3 
X\ + Y159 + Z\ + C/1345 

-^23 + ^26 + -^23 + ^12345 
X3 + I379 + Z2 + U24 
Xu + Y 48 + Z 13 + U b 

T,X + TY + + T.U 



5 — 



J 
k 
I 

m 
11 
-D/2 



(3.15) LAB: sistema-p2 



note that as the diagram does not have massive internal lines, the system have not "vari- 
ables" (sum indices) Wj. Indices Uj are concerned with Taylor expansion of, 



exp 



aj6 + adcj) + Bj6 + 



A 



+ ^9uj 2 

P2 



Such system has in principle 77,520 (Cj°) solutions. Determinant vanish in 42,526 of 
them, i.e., we must search hypergeometric type functions representing the original two- 
loop integral in a total of 34,994, less than half of possible solutions for the 7 x 20 system. 
However, in this case such series in far more complicated, since they are 13-fold ones. The 
results of such analysis will be presented elsewhere. 



4. Conclusion 



We studied in this work several integrals pertaining to non-planar double-box diagrams. 
Firstly we considered the case where four internal particles have mass fi, then write down 
the result in terms of a triple hypergeometric series. The second part deals with massless 
internal particles but off-shell external legs. We calculated the generating functional for 
our negative-dimensional integrals and then presented some particular cases of interest. 
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